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Abstract. Theoretical and computational advances in lattice calculations are reviewed, with focus 
on examples relevant to the unitarity triangle of the CKM matrix. Recent progress in semi-leptonic 
form factors for B nlv and B D*lv, as well as the parameter ^ in B^-B^ mixing, are highlighted. 



INTRODUCTION 



> 

^ I To test the CKM picture of CP and flavor violation, a combination of theory and 

' experiment is needed. A vivid way to summarize the need for redundant information 

is the unitarity triangle (UT), sketched in Fig. 1. It depicts two triangles, ByA, which can 
(in principle) be determined from tree processes (e.g., semi-leptonic decays), and aC|3, 
in which the amplitude for 5*^-5^ mixing is involved [1]. The "tree" triangle can, for the 
. sake of argument, be taken as a measurement of the CKM matrix. Then, the "mixing" 

P-i! triangle tests whether new physics must be invoked to explain B^-EP mixing. 

p I To obtain the sides A, B, and C of these triangles one must calculate hadronic proper- 

ly ■ ties from first principles of QCD. Here "calculate" implies that a reliable error estimate 

is given: one that includes all sources of uncertainty. Lattice gauge theory is well suited 
to the task: semi-leptonic form factors and the mixing matrix elements are conceptually 
^ ■ straightforward. Indeed, according to Martin Beneke [2], "[the] standard UT fit is now 

entirely in the hands of lattice QCD (up to, perhaps, iVif^,])." 

Until recently, lattice QCD has been burdened by something called the quenched ap- 
proximation (explained below). A bit of good news is that the available computer power 
is now sufficient to get rid of this approximation [3, 4, 5]. Another bit of good news is 
that several lattice groups have used the quenched approximation in the spirit of a blind 
analysis: although quenching could change the central value, one can analyze all other 
uncertainties as if the underlying numerical data were real QCD. This exercise has left 




FIGURE 1. Unitarity triangles: on the left is the "tree" triangle; on the right, the "mixing" triangle. 
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us with several quantities, including those needed in the UT fits, with full error analy- 
ses, apart from quenching. In the next few years, we should have full QCD calculations 
with a complete assessment of the uncertainties, and thereafter the uncertainties can be 
incrementally reduced. 

Because it is important to understand the theoretical uncertainties, this talk will start 
by sketching where they arise. Most reviews focus too much on central values, so, even 
when turning to numerical results, the focus here is on the error bars. 

The paper ends with a short swim through the "Octopus's Garden", that is, some 
recent work by Nathan Isgur, who passed a way a few weeks before Heavy Flavors 9 
convened. The symposium, and this paper, are dedicated to his memory. 

UNCERTAINTIES IN LATTICE CALCULATIONS 

Lattice QCD calculates matrix elements by integrating the functional integral, using a 
Monte Carlo with importance sampling. The Monte Carlo leads to (correlated) statistical 
error bars. This part of the method is well understood for quenched QCD and, these days, 
rarely leads to controversy. When conflicting results arise, they originate in different 
treatments of the systematics. The consumer probably does not need to know how the 
Monte Carlo works, but should develop an intuition of how the systematics work. 

The main tool for controlling systematics is effective field theory. In this talk, we 
are concerned mostly with three classes of effects: those connected with the lattice 
spacing a, the heavy quark mass(es) mg, and the light quark mass(es) m^. Inside the 
computer there is a hierarchy of scales 

niq <^ A<^ iuq, n/a, (1) 

although in practice 7l/a is a only several GeV, and the "light" quarks are never as light 
as the up and down quarks. The QCD scale A gauges the size of power corrections; ex- 
perience suggests it is 700 MeV, give or take a factor of -\/2. With familiar techniques of 
effective field theory, lattice theorists can control the extrapolation of artificial, numeri- 
cal data to the real world, provided the data start "near" enough. Similarly, non-experts 
usually have an intuition of how effective field theories work, so they can check, on the 
back of the envelope, whether systematic errors have been treated sensibly. 

The notable exception to the rule of effective field theory is the valence, or quenched, 
approximation. Consider the pictures in Fig. 2. The one on the left depicts a meson 
made of a valence quark and anti-quark, bound by a shmear of gluons. The gluons can 
also create virtual quark-antiquark pairs, leading to the picture on the right. These are 

-0 ^ -i^^ 

FIGURE 2. Quenched approximation: processes on the left are incorporated "exactly", whereas pro- 
cesses on the right are not computed, but modeled with a shift in the bare couplings. 
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computationally very costly. The quenched approximation omits them but compensates 
the omission with a shift in the bare couplings. This is analogous to a dielectric, where 
gQ ^q£. Quenching retains many effects, such as retarded gauge potentials, that are 
omitted in, say, the quark model. It is also the first term in a systematic expansion [6]. 

The quenched approximation can fall short of reproducing chiral logarithms of the 
form ln(A^/m^) [7, 8, 9, 10]. Figure 3 shows some quark-line configurations that 
generate, at the hadronic level, meson loops. The quenched approximation includes (b), 
but not (a) or (c). As a consequence, some pion loops are omitted, and r\' loops are 
mistreated. A better situation is a "partially quenched" theory, with dynamical quark 
loops, but, possibly, myaience 7^ 'Wsea- Then features like the ri'-7i splitting emerge, and it 
is possible to relate the partially quenched theory to full QCD using chiral perturbation 
theory. We shall return to chiral logs in the next section, when discussing B^-B^ mixing. 

To understand lattice spacing effects Symanzik suggested matching lattice gauge 
theory to continuum QCD [12]: 



The right-hand side is a local effective Lagrangian (LEz ) renormalized in some con- 
tinuum scheme; details of the scheme are not important. Discretization effects are, of 
course, short-distance effects, so as usual in an effective field theory, they are lumped 
into the coefficients Ki. 

There are two key points to the Symanzik effective field theory. First, if Ka is small 
enough, the can be treated as a perturbation. For example, for the proton mass 



using the leading term for Wilson fermions as an example. Here csw is the so-called 
clover coupling, and Kf^p = (1 — csw)/4 + 0(av). Second, to reduce lattice spacing 
effects, one can tune csw so that K^p vanishes or, in practice, is 0(af) or 0{a). Thus, 
the Symanzik effective field theory shows that lattice artifacts can be reduced through 
short-distance process-independent methods. Indeed, for light hadrons a combination of 
this Symanzik improvement and extrapolation in gives continuum QCD results with 
very well understood uncertainties (modulo quenching). 

With the bottom and charmed quarks, the mass is large in lattice units: m^a ~l-2 
and mchO- about a third of that. The split in the LEx between QCD and small correction 
breaks down, because there are terms in the in Eq. (3) that go like (niQa)". It will 
not be possible to reduce a enough to make m^a <C 1 for a long time: Moore's Law 



Aat = Xqcd + J^«'''i«^!(g^,m^a;//)0/(//). 



(2) 



mp{a) = mp - aKaF{csw){p\W ■ FMp) 



(3) 




(a) 



(b) 



(c) 



FIGURE 3. Quark line configurations that lead to meson loops. Adapted from Ref. [11]. 
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TABLE 1. Widely used methods for heavy quarks in lattice gauge theory. 



# method Ref. how HQET enters 



1. static approximation [13,14] 

2. lattice NRQCD [15,161 

3. extrapolation from below charm ad hoc 
3'. 3 + 1 ad hoc 

4. "Fermilab" [17, 18] 



-^^static = -hDA,h 

discretization of first few terms of HQET 
results with small niQ fit to Taylor series in l/niQ 
as in method 3 

match lattice to QCD, term-by-term in HQET 



suggests 15-25 years. There are, nevertheless, several ways to treat heavy-light hadrons 
in lattice calculations, all of which appeal to HQET in some way. Table 1 lists the most 
widely used methods. In lattice ]^QCD and in the Fermilab method, it is possible to set 
niQ = nih, even when ^ miy. There are, of course, uncertainties involved, but a grasp 
of the basics of heavy-quark theory is enough understand them. 

A convenient way to contrast the uncertainties in the various methods is to match 
lattice gauge theory to (continuum) HQET [18]. Instead of Eq. (3) one writes 

V- lat T, Tr. , hD^h hiH-Bh 

Lm = ^HQET = CnOn = -mihh - hD^h + + — + • • • . (4) 

The operators On on the right hand side are the same as in the HQET description of 
continuum QCD. The short-distance coefficients c]f', on the other hand, are not the 
same, because there are two short distances, a and 1/mg. Heavy-quark lattice artifacts 
are, thus, isolated into the mismatch c}^^ — In lattice NRQCD and in the Fermilab 
method, these mismatches can be reduced, along the same lines as reducing Kc^p in the 
Symanzik program. In recent work, they are controlled to several percent or less. 

IVIost work with method 3 has chosen normalization conditions for which the mis- 
match in the kinetic and chromomagnetic terms, although formally O(mga)^, is quite 
large in practice. It is also not well understood how this mismatch is amplified when 
extrapolating linearly or quadratically in 1/mg from mg ~ 1 GeV up to m^. 



LATTICE CALCULATIONS 

We now turn our attention to some of the most interesting recent lattice calculations in 
B physics. The discussion focuses on the error bars, and central values are deferred to 
the next section. We consider the ratio ^ for B-B mixing in unquenched QCD, and the 
semi-leptonic decays B %lv and B ^ D*lv for \Vub\ and \ Vcb\ in quenched QCD. 



Neutral Meson Mixing 

The mass difference of neutral B meson CP eigenstates is 



Am^o = ^^\V;M^B{B',\Qf='\B',) +new physics? (5) 
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where q G {d, s}, Sq is an Inami-Lim function, rig is a short-distance QCD correction, 



and ^ is a four-quark operator in the AB = 2 electroweak Hamiltonian. Note that 



'9 



new physics could compete with the Standard Model. The matrix element of is 
usually written 



AB-- 



8 



(6) 



but lattice QCD gives {B^q\Q^=^\B^^) and fe^^ (from (0|A^^|5j)) individually. This tra- 
ditional separation turns out to be useful when looking at chiral logs. 

Conventional wisdom says that the uncertainties in Bb^ and in = fB^BsJ/g^Bg^ 
are "easy" to control, because they are ratios. At Lattice '01, Norikazu Yamada of 
JLQCD reported on new results [19] that suggest otherwise. The JLQCD collaboration 
is mounting a large project to carry out calculation with the loops of nf = 2 flavors of 
quarks. Figure 4 compares their previous quenched work for /g^ with their new (and still 
preliminary) work with nf = 2.\n the spirit of a blind analysis, the definition of ^f^^ is 
not so important. The important matter is whether the curve is a straight line, or whether 
there is any curvature as a function of the light quark mass. (Because the pseudoscalar 
oc niq, the horizontal axis is, essentially, the light quark mass.) 

The quenched approximation (diamonds in Fig. 4) shows no evidence of curvature. 
With Hf = 2 (triangles in Fig. 4), however, allowing for curvature obtains a better fit. 
In fact, curvature is expected from chiral logarithms. Meson loops with fixed spectator 
mass Ms and varying sea quark mass m/ contribute to fg, [20] 



A/s, = co + cim 



11// 



l+3g 
(47rA) 



m: 



■ll.'/ 



(7) 



where rj^^/ denotes a pseudoscalar meson with constituents q and q' . From the D* width, 
~ 1.04 [21]. The second plot in Fig. 4 shows that Eq. (7) describes the points well. 
It is worth making a few remarks about this kind of analysis. The formula for /g^ 
(where the spectator and sea quarks have the same mass) is slightly different, so it is 



M ~ 6 GeV(~B meson) 




0.2 0.4 0.6 



2 2 

m [GeV ] 



PQyPTfit with fixed log coef. 



heaviest ^ 



• B meson 



X/dof=0.32 
X/dof=0.12 



3 4 5 

,2 



FIGURE 4. Recent calculations of the decay constant 4>b = fB^/mB from JLQCD. 



9th International Symposium on Heavy Flavor Physics 



September 10-13, 2001 



5 



difficult to display all chiral log effects on one plot. Nevertheless, one can see from 
Fig. 4 that the effect is to increase fsjlsj- The chiral logs for Bg^ are multiplied by 
1 — 2)g^ pa 0.04, so the BbJBb^ is almost insensitive to them. Therefore, the chiral logs 
in the ratio ^ come almost completely from faJfEd- One concludes that ^ may have 
been underestimated in the quenched approximation. A new estimate, based on JLQCD 
and previous work, is given below. 



B nlv and Vut 

The semi-leptonic decay B — > nlv is mediated by a — > w transition. The decay rate is 



dE 127i3 



-%i,\'\U{E)\\ (8) 



where v = ps/mB is the B meson's velocity, E = v- pji, and = E^ — m^. The form 
factor f+{E) is a linear combination of the form factors f±{E) and f\\{E), defined 
through the matrix element of the vector current 



{%\V^'\B) = yf\\{E) +pIME)] . (9) 

The pion energy E is related to = + m^. — InisE. Chiral symmetry and heavy- 
quark symmetry are simpler to follow with /|| and f± considered as functions of E, 
rather than /+ considered as a function of q^. For example, heavy-quark symmetry 
suggests relations between B and D decays with the same E. In principle, one would like 
to compare the E dependence of experimental measurements with lattice calculations. 
Discretization uncertainties are smallest for low E, where phase space suppresses the 
event rate. Therefore, lattice calculations and experimental measurements will have to 
be combined in the way that minimizes the error on \Vuh\. 

In the past year or so, four new quenched calculations of /+ (E) for B nlv appeared, 
using several different methods [22, 23, 24, 25]. It is, thus, timely to compare and con- 
trast. In addition to discretization effects at large pa, there is evidence for considerable 
dependence on the light spectator quark mass [24]. 

The calculations of UKQCD [22] and APE [23] use method 3, so they have 
discretization errors of order (mga)^. They 

keep mg < 1.2mc and extrapolate up to nib TABLE 2. Heavy quark masses used in 
linearly or quadratically in 1/mg. The quark Ref. [22]. Note that mp^jg < m^. 
masses used in Ref. [22] are shown in Table 2. 
Those in Ref. [23] are a bit larger, which is 
better for HQET but worse for the discretiza- 
tion errors. In these papers models are used to 
extrapolate to the full kinematic range of £. It 
is better to think of these calculations as com- 
puting the model parameters in D decays, and 
then invoking HQET to make predictions for 



K 


aniQ 




Mp 






(GeV) 


(GeV) 


0.1200 


0.485 


1.52 


2.035(5) 


0.1233 


0.374 


1.23 


1.771(5) 


0.1266 


0.268 


1.02 


1.483(5) 


0.1299 


0.168 


0.69 


1.157(5) 
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B decays. The error associated with the assumptions are, at least to me, not transparent. 

The calculation of El-Khadra et al. [24] uses the Fermilab method, and the calculation 
of JLQCD [25] uses lattice NRQCD. Since the matching procedures build in the heavy- 
quark expansion, both calculate directly at mg = m/,. To avoid introducing models of the 
E dependence, these two papers advocate comparing lattice and experiment in a region 
where the two overlap. For example, one can look at 

fEmax ^ ^11 ?TT^ fEmax rIT 

Tb{E^..) = / dEp'\ME)\^ = Wl27^ / ^' (10) 

Jm^ \Vubr GpniB Jm^ dE 

where ^max is an upper kinematic cut. At present ^max ~ 1 GeV, but, with future 
increases in computer power, it may be possible to raise the cut. This method is sketched 
in Fig. 5 using the form factor from Ref. [24]. Uncertainties are still about 15-20% 
on \ Vub\- A strategy for reducing them are given in Ref. [24]. 

Ref. [24] finds strong dependence on the light quark mass, as illustrated for 
f±{OJ GeV) in Fig. 5. This is the only calculation to reduce the light quark down 
to niq — OAm-s. The extrapolation to physical quark mass is the largest source of uncer- 
tainty quoted by Ref. [24], and a similar uncertainty is presumably present in all four 
calculations. One would like even lighter quarks, which take more CPU time. 

Because of the discretization errors at higher E, the large light quark mass depen- 
dence, etc., a sanity check on f+{E) for B nlv would be welcome. The calculation of 
similar form factors in D decays encounters many of the same issues. A compelling test 
would be to compare lattice and experiment for /f^^///), and /f ^"'//d, as a function 
of E. New physics is unlikely to alter the underlying processes, and the CKM matrix 
drops out. The CLEO-c program [26] promises to measure these ratios with an uncer- 
tainty of a few percent. 




FIGURE 5. B ^ nlv: (a) Integrand of Eq. (10). (b) Dependence of the form factor /x on the light quark 
mass, in the range QAnis < nig < 1.2ms . 



9th International Symposium on Heavy Flavor Physics 



September 10-13, 2001 



7 



B D*/v and Vcb 

The exclusive semi-leptonic decay — > D*/v is a good way to determine \ Vcb\, but an 
estimate of the hadronic transition is needed. The differential decay rate is 



dT 

dw 



47C3 



^ 

^ V - Iml,* (niB - mo* f- Q {w)\Vcb\^\ !Fb-^d* (w) 



(11) 



where w — v-v' and ^ (1) = 1. At zero recoil (w — I) heavy-quark flavor symmetry 
forbids terms of order 1/mg. The exclusive determination of \Vcb\ therefore proceeds 
by extracting \ Vcb!FB-^D*{.^) \ from a fit to the measurement of dT/dw, and then taking a 
theoretical estimate of jfB^£)*(l). One would prefer the latter not to depend on models. 

To see what is needed from the theoretical calculation, let us review the anatomy of 
!Fb^d*{^)- From HQET one can show that, at zero recoil. 



!Fb^D* ( 1 ) = T|a llsgur-Wise + Omke/ ^ + 5i/^2 + 5 



(12) 



where r[A is a short-distance radiative correction, and the 5i/„n contain the long-distance 
properties of the bound states. The numerical values of r[A and the 5i/„« depend on how 
HQET is renormalized. It is less important which renormalization scheme is chosen than 
it is to use the same scheme for both. The 1 /m^ corrections take the form 



Sl/m2 — 
Sl/m3 = 



+ 



2£a 



{InicY {lmc){lmb) {Inib) 



(3) 



+ 



(3) 



(2mc)3 (lmc){2mb) {ImbY 



(13) 
(14) 



where E = 1/ {Inic) -f- 1/ (2m/,) and A = 1/ {Inic) — 1/ {Inib). 

With brute force alone, a sufficiently precise calculation of U^b^d*{^) lies beyond 
reach. (See Ref. [27] for details.) Hashimoto et al. [27] have devised a method to extract 

all £sinEq. (13) and all but £^3'^ inEq. (14). The key is the heavy-quark mass dependence 
of three double ratios of matrix elements. As in Eqs. (12)-(14), HQET implies 



{D\ci^b\B) {B\bfc\D) 
{D\cfc\D){B\bfb\B) 
{D*\cfb\B*){B*\bfc\D*) 
{D*\cfc\D*){B*\hfb\B*) 

{D*\cyjy5b\B){B*\byjy5c\D) ^ 

{D*\cyJy5c\D) {B*\byJy5b\B) r]A^cr]^bb 



^Acb^Ahc 



I -A' 



(15) 



(16) 



(17) 



In particular, heavy-quark spin symmetry requires the same is to appear in Eqs. (15)- 
(17) as in Eqs. (13) and (14). 

In a lattice calculation of these double ratios most of the statistical and sytematic un- 
certainties cancel. The main difference from continuum QCD is that the short-distance 
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coefficients are different [18, 27]. But, just like in continuum QCD, the short-distance 
behavior can be computed in perturbation theory. Thus, the analysis proceeds by remov- 
ing the lattice short-distance contribution, fitting to the prediction of HQET, and then 
reconstituting (1) from the £s and r\A- The result is 

cr (U—ri QO9+0.024 I A ni ^;+0.003+0.000+0.006 on 
iFfi^DHl) -U-y-' •'^-0.017 ±^-^l"-0.014-0.016-0.014- (1^) 

To encourage the reader to think about the uncertainties, the central value is not revealed 
until the end. The error bars come from, in order, statistics and fitting, matching lattice 
gauge theory with HQET to QCD, lattice spacing dependence, light quark mass effects, 
and the quenched approximation. The secret to the small error on f is that they all scale 
as ^ — 1, by design. As a fraction of J — 1 the uncertainties are still sizable: 5-25%. 

The most novel aspect of this analysis is how seriously it takes the idea of matching 
lattice gauge theory to QCD through HQET. A central part of the analysis is to calculate 
short-distance properties, which Ref. [27] does partly in perturbation theory. The sec- 
ond error bar reflects the associated uncertainty. It is reducible, but through traditional 
theoretical physics, rather than intensive computation. 



AN OCTOPUS'S GARDEN 

Although it does not bear directly the unitarity triangle, I would like to discuss some 
recent work by Nathan Isgur. I did not know Nathan well, but his enthusiasm for physics, 
especially the strong interactions, made a big impression on me. He struck me as the kind 
of person who had a hand all sorts of things. He must have needed eight arms to keep all 
his projects moving. 

In the last year of his life, one of his projects was to understand whether instanton- 
like gauge fields or, on the other hand, disordered gauge fields play a larger role 

in QCD [28, 29, 1 1]. This work reflected Nathan's broad knowledge of the phenomeno- 
logical and theoretical sides of the strong interactions, touching on the OZI rule, chiral 
zero modes, the r]' mass, instantons, the quark model, the AdS/CFT correspondence, 
lattice QCD, and the large Nc limit. A particularly striking passage notes first how recent 
work by Witten on the AdS/CFT correspondance favors disordered, confining gauge 
fields (although those sentences were probably written by Nathan's collaborator Hank 
Thacker), and then how details of how to treat strong decays in the quark model favor 
the disorder also. That part was certainly written by Nathan. 

These issues surround the quantum ground state, or vacuum, of QCD. Gauge theories 
have many classical ground states, one for each integer. Instantons are the semi-classical 
configurations that tunnel from one ground state to another. In a quantum mechanical 
situation, fluctuation-dominated configurations could also mediate tunneling. Nathan 
proceeded by devising tests that could distinguish whether the quantum vacuum is 
obtained via such disordered gauge fields or via a gas or liquid of instantons. He then 
used (quenched) lattice QCD to see which way the gauge fields behaved. 

One test stemmed from the observation, from empirically successful quark models of 
strong decays, that quark pairs pop out of the vacuum with scalar (i.e., ^Pq) quantum 
numbers. That implies that the OZI rule should fail for 7^^ = 0++ quantum numbers. 
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(The OZI rule says that decays, like (|) nn, in which quarks annihilate, are weaker than 
those, like (|) — > KK, in which they do not.) Empirically, the OZI rule fails in the ^ 
sector; that is how the T|' mass is split away from the other pseudoscalar mesons. For the 
0"^"*" sector, however, there is a dearth of experimental data. Isgur and Thacker's idea to 
to compute the quark annihilation process with lattice QCD, for various J^'^, which is 
possible at leading order in the OZI approximation even in quenched QCD [28]. They 
found a small amplitude where the OZI rule holds empirically, e.g., for J^^ = 1~+. But 
for scalar and pseudoscalar channels the amplitude is large. This is evidence against 
an instanton-dominated vacuum, because instantons couple to quarks through the axial 
anomaly, that is, preferentially to the pseudoscalar channel. 



UNBLINDING 

Instead of a paragraph of bland conclusions, the reward for readers who have made it this 
far consists of the most interesting numerical results. The preliminary results of JLQCD 
on B physics with n/ = 2 are [19] 

fB, = 190(14)(7)(19)MeV, (19) 
IbJIb, = 1.184(26)(20)(15), (20) 
i = 1.183(27)(20)(15), (21) 

where ^ is a combination of their results for /bJ/bj and BbJEb^.I^o. her review at Lat- 
tice '01 [30], Sinead Ryan took stock of the ratio ^, which enters into UT fits. Although 
JLQCD 's result is still preliminary, it cannot be denied that the chiral logarithms could 
affect fsJfBd hence, ^. Ryan's average, which strikes me as reasonable, is 

^=1.15±0.06lj;:j!J (22) 

which retains the central value in common usage, but allows for a future upward revision, 
once the chiral behavior is fully understood and controlled. 

For the zero-recoil form factor needed to determine \Vcb\ from B ^ D*lv, Ref. [27] 
finds 

J,_z)Kl)=0.913lSiS, (23) 

with systematics added in quadrature. Note that this result does not include the QED 
correction of +0.007, which is included, for example, in the BaBar Physics Book. This 
result still relies on the quenched approximation. Nevertheless, it is probably still under 
better control than estimates based on the quark model. 

A hallmark of these two results, shared by some work on 5 ^ nlv for \Vub\ [24, 25], 
is that they attempt a full analysis of the uncertainties, including those of heavy quarks, 
the chiral behavior, and quenching. Thus, they are suitable templates for the next round 
of lattice calculations, from which one can expect serious unquenched calculations with 
a direct impact on our knowledge of the unitarity triangle. 
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